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Abstract 

We propose a consistent quantization of the Berry-Keating Hamiltonian xp, which is currently 
discussed in connection with the non trivial zeros of the Riemann zeta function. The smooth 
part of the Riemann counting formula of the zeros is reproduced exactly. The zeros appear, not 
as eigenstates, but as missing states in the spectrum, in agreement with Connes adelic approach 
to the Riemann hypothesis. The model is exactly solvable and renormalizable, with a cyclic 
Renormalization Group. These results are obtained by mapping the Berry-Keating model into the 
Russian doll model of superconductivity. Finally, we propose a generalization of these models in 
an attempt to explain the oscillatory part of the Riemann's formula. 

PACS numbers: 02.10.De, 05.45.Mt, ll.10.Hi 
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I. INTRODUCTION 



The Riemann hypothesis (RH) is a central problem in Pure Mathematics due to its 
connection with Number theory and other branches of Mathematics and Physics. The RH 
is the statement that all the non trivial zeros of the zeta function C(s) lie on the critical line 

n 

Re(s) = 1/2 .1]. Most of the physical approaches to prove the RH are inspired by the old 
Polya and Hilbert conjecture, which states that the imaginary part, E a , of the Riemann zeros 
~ + iE a ) = 0, are the eigenvalues of a Hamiltonian and thus real numbers if the RH is true 
. This approach is supported by the statistical properties of the zeros, the Montgomery- 



Odlyzko law 



based on the Gaussian Unitary Ensemble distribution (GUE) [a], and 



also by the quantum chaos interpretation of the oscillatory part of the Riemann counting 
formula of the zeros j^j]. Berry has proposed that the Riemann dynamics is given by a 
classical chaotic Hamiltonian with isolated period orbits, whose quantization would yield a 
point like spectrum containing the E' a s jf]]. Other physical approaches are based on statistical 
mechanics ^sj], superconformal invariance j^], supersymmetric quantum mechanics lfl| . etc 
sec . Il2j ] for recent reviews). An interesting related approach is the construction of a 
quantum mechanical potential containing in its spectrum the prime numbers jisl ]. 

The starting point of our work is the Hamiltonian i^BK — %P, proposed by Berry and 
Keating, which reproduces at the semiclassical level the smooth part of the Riemann's 
formula giving the number of non trivial zeros below a given number, E QQ. ths 
Hamiltonian is formal since a consistent quantization has not yet been found. In this paper 
we shall propose a solution of this problem defining xp on a lattice. This is achieved working, 
not with xp, but with its inverse 1/xp, which turns out to be related to a class of QM models 
with limit cycles, or rather centers, in the Renormalization Group (RG). 

The idea that the RG may have limit cycles was first considered by Wilson in 1971 
however at the time no models with this behavior were known. In the last few years limit cy- 
cles in the RG have been discovered in various models in several physical contexts, including 

" lsl llflj ] , quantum mechanics |2I 



nuc 
ity 



ear physics 



21 



13], quantum field theon 



, superconductiv- 

22|, Bose-Einstein condensation [23], effective low energy QCD |24| . S-matrix models 



25|,|26(, few body systems and Efimov states [17|, |27|, etc (for a review of see [27[). The sub 



ject of duality cascades in supersymmetric gauge theory 28] is also suggestive of limit-cycle 
behavior. The possibility of chaotic flows has also been recently considered fsoL l^. 
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The relevant model for this work is the, so called, Russian doll BCS model of super 



conductivity, and specially its QM version |21l. I22j. In the RD model the standard pairing 
coupling g flows periodically with the scale, a fact that is intimately related to the exis- 
tence of a series of bound states of two electrons (Cooper pairs) whose energies scale as 
e -nA ^ n _ o ; 1, 2, . . . ), where A is the period of the RG cycles. These bound states have a 
size that scales as e nX (n = 0, 1, 2, . . . ) which is the reason for calling them Russian dolls, by 
analogy with the popular matrioskas. The RG period is the main feature of a RG with limit 
cycles. In the RD model it is given by A = 2ir/h, where h is a coupling in the Hamiltonian 
that breaks the time reversal symmetry. 

Why should this model be related to the Hamiltonian xp which seems so far apart? Some 
hints lie on the following observations. In the RD model the wave function of a Cooper 
pair with energy E is given approximately by %p{n) ~ l/(n — E) l ~ lh , where n — 1, 2, . . . , N 
label the energy levels of the electron pairs. Since in this model the energies E n converge 
towards zero, E n ~ e~ nX , there exist a bound state with E = 0, whose wave function is 
i/j(n) ~ This form recalls the Dirichlet series of the zeta function, ((s) = l/n s , 

with s = 1 — ih. The latter sum is the interaction term that is multiplied by the coupling g. 
Apparently, the RD model should be connected to ((1 — ih) rather than to <^(l/2 — iE). It 
is worth to mention that in the cyclic sine Gordon model, which also has a cyclic RG, the 
zeta function £(1 — ih) appears in the expression of the effective central charge, where h is 
related to the period of the RG cycles as A = n/h jlfll. I^. 

The situation is different in the BK model where the formal eigenfunctions of the normal 
ordered Hamiltonian (xp+px)/2 are given by ip(x) = l/x l l 2 ~ tE , that resembles the Dirichlet 
series of ((1/2 — iE). Compare to this, the RD wave function with E = has in the exponent 
of n the factor 1 instead of 1/2, and the imaginary part is fixed to a constant h instead of 
being energy dependent. We shall show in this paper that these two problems can be solved 
at once by relating the BK and the RD models by means of a third model whose classical 
Hamiltonian is simply the inverse of the BK one, namely 1/xp. In a deep sense, the BK and 
the RD models turn out to be the two faces of the same coin. 

The organization of the paper is as follows. In section II we review the BK and the RD 
models. We also define the inverse model based on 1/xp and discuss its relation with the 
previous ones. In section III we solve the inverse model in the continuum limit. In section 
IV we study its renormalization showing the existence of RG cycles and their connection 
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with the continuous solution. In section V we give its exact solution, derive the smooth 
part of the Riemann counting formula and present a numerical analysis. In section VI we 
generalize the inverse model and begin to study its properties. We state the conclusions and 
prospects in section VII. 

II. THE HAMILTON! ANS 

In this section we introduce the Hamiltonians of the Berry-Keating model, the inverse 
model (X) and the Russian doll model, and study their relations which can be represented 
symbolically as 

BK^l=BK~ l ^RD (1) 
Before presenting the details, we shall start with an overview. 

As we explained above, the classical version of the Berry-Keating Hamiltonian Hbk is the 
product xp, where x and p are the position and momentum of a particle moving in one 
dimension. Hbk fullfills some of the asumptions of the quantum chaos approach to the RH, 
particularly the breaking of time-reversal symmetry, to accomodate the GUE hypothesis 
and suggestive analogies involving the zeta function, the trace formula, the Riemann-Siegel 
formula, etc. The BK proposal has however remained at a speculative level since a consistent 
QM model of xp has not been constructed in so far. 

To solve this problem we shall define the Hamiltonian Hj, which at classical level is the 
inverse of the BK Hamiltonian, namely 1/xp. It would seem that nothing is gained by this 
trick. Nevertheless, a consistent QM model can be constructed by quantizing 1/xp on a 
lattice. The key observation is that the inverse of the momentum operator p = —ih d/dx is 
the ID Green function j^G(x, x') = ^sign(x—x'), where sign(:r) is the sign function. We shall 
define a regularized lattice version of the operator p" 1 as the matrix = ^ sign(n — m), 
where n and m run over the integers 1, 2, . . . , N . The quantization of the position operator 
x is the diagonal matrix X n m = n<5„ im . In this construction the Hamiltonian Hj acts on a 
discrete Hilbert space of dimension N. 

It will be important to study the behaviour of Hj under the Renormalization Group 
transformation which integrates the highest energy level reducing the size of the system to 
N — 1. For the model to be renormalizable one needs an extra term with coupling constant 
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go The other coupling constant, hi = 1/h, multiplies the term containing Hence 
the X model depends on two couplings constants hi and gi. Under the RG transformations 
hi remains invariant, while gi flows periodically with the scale. This periodicity is in turn 
related to the spectrum of the X model. 

The third Hamiltonian is the QM Russian doll Hamiltonian, which is a generalization 
of the standard BCS model of superconductivity having a cyclic RG. The RD Hamiltonian 
Hud also depends on two coupling constants go and ho, where ho multiplies a time- reversal 
breaking term, proportional to P~m, while go multiplies the familiar s-wave pairing inter- 
action of the BCS model which preserves the time-reversal symmetry. The spectrum of this 
model has a series of bound and antibound states with Russian doll scaling. The RD model 
is brought in due to its close relationship with the X model. Indeed, we shall show that each 
eigenstate of the X model coincides with the zero energy bound state of an associated RD 
model. This bound state appears at the threshold of that model. From the RG viewpoint 
these two models will also be related. 

A. The Berry-Keating model 



The classical BK Hamiltonian 



14, 



Hi K = xp, (2) 
has classical trayectories given by the hyperbolas 

x(t) = x e\ pit) = p e~*. (3) 

The dynamics is unbounded and one should not expect a discrete spectrum at the quantum 
level. Despite of that, Berry and Keating regularized the model introducing the Planck cell 
in phase space with sides l x , l p and area h = l x l p , such that \x\ > l x and \p\ > l p They 
computed the number of states Nbk{E), with an energy below E, using the semiclassical 
formula 

A(F) 1 r E i lp r E i x 

N BK (E) = ^ = - dx dp+... (4) 

J l x J l v 
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where A(E) is the area of the region l x < x < E/l p , l p < p < E/l x , xp < E (see fig. la). 
The result is 

N BK (E) = ^ (log + (5) 

Taking into account a Maslov phase —1/8, BK finally obtain 

where the energy E is measured in units of h. Eq.© agrees with the asymptotic expansion 
of the smooth part of the Riemann counting formula of the non trivial zeros of the zeta 
function 

N sm (E) = hmlo g r(]+ i -E) --^logTT+1. (7) 

7T y4 A J Z7T 

This expression is tt" 1 times the phase of the zeta function £(1/2 + [30|. The non trivial 
zeros, p, are known to lie in the critical strip < Re(p) < 1. The number of them, N(E), 
in the range < Im(p) < E, is given by the Riemann formula, 

N(E) = N sm (E) + N osc (E) J (8) 
where the oscillatory part depends of the zeta function on the critical line, 

N OBC (E) = -lm\og((l + iE) (9) 



7T °'V2 

and it is of order \ogE j^jj. The derivation of eq.(JBJ is heuristic and, to our knowledge, 
has not been reproduced from the quantization of the BK Hamiltonian. Bhaduri et al. [20( 
obtained eq.((7j), up to some factor, studying the inverted harmonic oscillator Hamiltonian 
p 2 — x 2 , which is related by a canonical transformation to xp. They employed a phase shift 
approach instead of the semiclassical bound state approach. 

The BK work was motivated by Connes's work to prove the RH using p-adic numbers. 
This is known as the adelic approach to the RH and it is based on the construction of an 
abstract space where acts an hermitean operator, whose eigenvalues are the non trivial zeros 
of the zeta function The truth of the RH lies in the proof of a certain classical trace 
formula. Connes has also considered the operator xp in the adelic theory. The regularization 
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FIG. 1: la) The region in shadow is the phase space considered by Berry and Keating in the 
semiclassical computation of the number of states of the xp Hamiltonian (see eq.©). lb) The 
shadow region gives the phase space considered by Connes (see ea. (fTT)|0 . 

is not the Planckian cell, but a standard cutoff in phase space given by \x\ < A, \p\ < A. 
The number of states in the region < x < A, < p < A,xp < E (see fig. lb) 



2E - [ — + / dx dp + 

A/ Je/a Je/a 



(10) 



yields, 



iV Co ( £ ) = |lo g ^-|(l„ g |-l) + ... (11) 

where E is measured in units of h. There are important differences between eqs. ljllj) and Q. 
In Connes's formula the number of states diverges with the cutoff. There is also a negative 
sign of the term in pip, associated to the asymptotic Riemann formula, as compared with 
the BK formula (JBj) . According to Connes, this negative sign agrees with the overall sign in a 
brmal expression of N osc (E) obtained using the Euler product formula of ((1/2 + iE) into (jHJ) 
3ll | . These results imply that the Riemann zeros, in the adelic approach, are missing states 
in a continuum spectrum. The physical picture is that of white light with dark absortion 
lines labelled by the zeros of the zeta function. Using this analogy, the Riemann zeros, 
according to BK, would be emission lines. Our results are along the lines of Connes, as we 
shall show. 
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We shall end this brief review of the BK approach by giving the formal hermitian operator 
which corresponds to (J2J), 



H B k = \{xp+px) = -iH + J ) . (12) 



2 y r r ' \ dx 2 

whose formal eigenfunctions are 

M*) = xl/2 _ iE/n - (13) 

As noticed by BK, the wave function contains the power x~ s which appears in the Dirichlet 
series for ((s) and the Euler product, although those formulas are convergent only in the 
region Re(s) > 1, while in eq. (jlH)l s = 1/2 — iE/H, which lies outside. The rest of the BK 
paper discusses possible quantum boundary conditions that could generate the Riemann 
zeros. The guiding idea is that xp is the generator of the scale transformations in phase 
space, a fact that also plays a role in Connes's work. 
Generalization of the BK model 

The Hamiltonian (j2J) can be generalized to the form v(x)p, where v(x) is a generic function 
which we shall assume is positive and monotically increasing for x > 0. The classical evolu- 
tion equation x = v(x), implies that v(x) is the velocity of the particle. The semiclassical 
formula for the number of states a la Connes gives 

Nco(E) = ^ T + ±v-\E/A), (14) 

2vr J v -i iE /A) v{x) 2vr 

where t> -1 is the inverse function of v(x). If v(x) = x, eq.(|14|) becomes eq.pH). Similarly, 
the formal Hamiltonian associated to vp 

H B k = 7^{v{x) p + pv(x)) = -ih (^ v ( x )-^ + \ V '( X )J ' ( 15 ) 
has formal eigenfunctions 

. , v A (iE [ x dx' \ 

m ' ) = ms^[tLmP))- (16) 

If v(x) = x we recover JT 



S 



B. The inverse Hamiltonian 



The inverse of the generalized BK Hamiltonian (|15|) is 

HbIc = ^ 1/2 P- 1 (17) 
where we use the normal ordering v 1 ^ 2 p v 1 ^ 2 of Hbk, which is equivalent to (|15p. By 
asumption v (x) > for x > 0, then the square root of v(x) is real. For this and other 
reasons will shall work on the half line x > 0. The inverse of the momentum operator p is 

1 1 

{x\p~ l \x') = —G(x,x') = — sign(x — x'), (18) 

where G(x, x') is the ID Green function associated to d/dx and sign(x) is the sign function. 
The advantage of using instead of Hbk is that the former operator can be easily 

regularized on a ID lattice due to the simplicity of the ID Green function. The lattice 
version of the formal Hamiltonian ()17|) is the iV-dimensional matrix 

(n\Hi\m) = -f n (gi + i/i 7 sign(n - m)) f m , (19) 
where n, m — 1, 2 . . . , N and 

hi = p /„ = (20) 

Hj contains an extra term with coupling constant gj, which is needed for renormalization 
(see section IV). The pure BK model corresponds to the case gi = 0, but we shall also 
consider non vanishing values of gi. The regularized BK Hamiltonian will then be defined 
as the inverse of Hi, namely 

H^\v n , gi ,h) = Hr\f n ,gi,hi) (21) 
In the case where gi = 0, one gets 

H^\v n , gi = 0,h)= v{Xf' 2 P v(Xy/ 2 , (22) 
where X and are the N dimensional matrices 



(n\X\m) = n5 n>m , (n\P~ 1 \m) = — sign(n — m). (23) 
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FIG. 2: Eigenvalues of the commutator [X,P]/i (h = 1) for N = 50. We display 40 eigenvalues 
which converge towards 1. The remaining ones are larger in absolute value. The sum of all the 
eigenvalues is zero as a consequence of the vanishing of the trace of [X, P] . 

For N an even number, the matrix P~ l is not singular and its inverse is given by 



X and P do not satisfy the canonical commutation relation [x,p] = ih. The trace of the 
LHS of the commutator is zero while on the RHS is ihN. This is the textbook argument to 
show that the Heisenberg's indetermination relation cannot be realized by finite dimensional 
matrices. Nevertheless, most of the eigenvalues of [X,P]/i converge towards h (see figEJ) - 
From this result we expect to recover the formal BK model in the continuum limit of the 
discrete Hamiltonian (|22|). Since [X, P] ^ ih, the Hamiltonian ((221) cannot be written as 
l/2(v(X)P + Pv(X)), but this fact is unimportant. We can then work either with the 
Hamiltonians ()19j) or (j22j) - We shall choose the former for convenience. 

The gi interaction has a classical version. This term is non local in x-space and ultralocal 
in p-space. By dimensional reasons it should be proportional to S(p), yielding Hf = (l/p + 
5(p)gi/2)/v(x), whose inverse is again pv(x) (recall that p 5(p) = 0). At the quantum level 
gi does play a role in H BK (gi) = Hj l (gi) setting the boundary conditions. 

The operator K = h/xp was used by BK to implement a canonical transformation x — > 
x\ = h/p, p — > pi = xp 2 /h, named quantum exchange by the h dependence jl4| . Berry and 
Keating, tried to combine the dilatation and the quantum exchange symmetries, to generate 




(24) 
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FIG. 3: Pictorial representation of the RD model. The horizontal lines represent the energy levels 
£n ( n = 1, ■ • • , N). The dots are pairs of electrons occupying a doubly degenerate state. The arrows 
give the transitions among levels induced by the interactions in the Hamiltonian (|25[). 

the Riemann zeros. In our case, the operators xp and 1/xp are used as Hamiltonians and 
not as symmetries. 

Hi has two properties which simplifies its study: it is renormalizable and exactly solvable 
for all values of the coupling constants hj, gi and f n . One can find an analytic form of 
the exact eigenvectors and an explicit equation for the eigenenergies. These properties are 
reminiscent of the Bethe ansatz, as it is indeed the case. 



C. The Russian doll Hamiltonian 



A Hamiltonian closely related to Hi is 21] 



(n\H RD \m) 



- (g D + ih D sign(n-m)), 



(25) 



where £„,<?£>, hp are real parameters. In physical applications e n , are the kinetic energies of 
pairs of electrons occupying doubly degenerate levels in the conduction band of a metal. The 
term proportional to gu is a pairing interaction due to the phonon exchange in the s-wave 
channel, which leads to the formation of Cooper pairs if go > 0. The model with hp = 
was first considered by Cooper and it was the precursor of the Bardeen-Cooper-Schrieffer 
model of superconductivity |2i 



The coupling hp was introduced in reference 



2l| . motivated by the work of Glazek and 



Wilson, who proposed an extension of the well known 2D delta function potential in order to 
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show the existence of RG limit cycles in a simple QM model [20]. The QM Russian doll model 



(|23jl. together with its 



detail in references 



many body version, also have a cyclic RG, which has been studied in 



21 



22|. The exact solution was obtained in |3J] using algebraic Bethe 
ansatz methods, showing that the RD model is nothing but an inhomogenous XXX vertex 
model with a boundary operator. In the vertex formulation the coupling hp parameterizes 
the quantum Yang-Baxter matrix, the coupling go parameterizes the boundary operator and 
the energy levels e n give the inhomogenities. The most natural choice of the energy levels is 
e n = n, in units of twice the electronic level spacing. The equally spaced BCS model for a 



finite and small number of energy 
of ultrasmall metallic grains (see 



evels have been studied intensively due to the fabrication 
~l 

35| for a review). This model is also relevant in Nuclear 
Physics, where it is known as the picket fence model, and in other potential applications 
in Quantum Optics, and dilute Fermi-Bose gases (see j3] for a review). Up to date it is 
not known wether the RD-BCS model of superconductivity explains some physical system 
existing in Nature or in the Lab. A possible reason is that the Hd coupling breaks the time 
reversal symmetry, while the microscopic laws of Nature do not. However, this symmetry 
can be broken spontaneously or explicitely by the action, for example, of magnetic fields. 

We shall next show the relation between the RD and the I models. Let's first compare 
their Schrodinger equations, 

H I tP = E I V, H RD <p = E RD <f), (26) 

which read explicitely 

Eiipn = \ Em=i fn(9i + ihi sign(n - m))f m ij; m (27) 
(e„ - E RD )(j) n = \ J2Z=i(9d + ih D sign(n - m))(j) m . (28) 

Given an eigenstate, tp, of the X model let us define the wave function 

0n = Mn, (29) 



writting (|2*7)l as 



I 1 N 

-^Ei<\) n = 7^^2(91 + ihi sign(n - m))0 m . (30) 



^2 th 2 
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If Ei 7^ 0, and dividing eq.()30p by Ei, one obtains eq.()28p with Erd = 0. This leads to the 
identifications: 



£ n = 72' 9d = ~et, h D = — , (31) 

In h I h I 

which establish a map from the spectrum of the 1 model into a collection of RD models 
which have in common the energy values, e n , focusing on the Erd = state, while varying 
the couplings go and ho according to Ej. Moreover, eqs.(J20J) establish a link between the 
regularized BK model and the RD model, Combining (|3TJl with (J2UJ) we get 

e n = v n , h D = — , (32) 

where E — 1/Ej is an eigenvalue of Hbk = Hj 1 . Let us remark that the transformation (f2T)|) 
is non unitary. For the model e n = n and in the limit N — > oo, <p n will be a normalizable wave 
function, corresponding to a bound state, but tp n will not, as corresponds to a scattering 
state. 

The map BK — > RD is quite remarkable. It implies that the spectrum of the BK Hamil- 
tonian can be found by looking at the zero energy bound states of a RD model where the hp 
coupling is fine tuned according to the energy of the state. In this mapping, the velocities 
v n become energy levels e n . In particular, the choice v = x corresponds to the equally space 
model e n = n, which is the common choice for physical applications 0Q- 

The correspondence (J32|) can be easily derived at the classical level. Take go = for 
simplicity. The classical energy (|23|) is E RD = e{x) — h D /p (with h = 1). Setting E RD = 
one finds p e(x) = hp, which is the classical energy of the BK model with v (x) = e(x) and 
E = h D . 

In summary, we have constructed a consistent quantum model which generalizes the BK 
Hamiltonian xp by means of the inverse model /, whose spectrum can be mapped into the 
zero eigenstates of associated RD models. As we shall see this connection is the key of the 
renormalizability and solvability of the BK model so defined. 

III. CONTINUUM LIMIT 

The three models described in the previous section can be solved exactly. However it is 
worth to solve them first in the limit where the lattice size iV is very large. In this limit the 
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discrete variable n will be considered as continuous and varying in the interval 1 < n < N . 
The Schrodinger eq.(j2Zj) becomes, 



N 



E iip {n) — — I dm f(n)(gi + ihi sign(n — m))f(m)ijj(m), (33) 
2 Ji 

where ip{n) and f(n) are continuous functions. To solve we use again the change of 
variables (J2HI), 

1 f N 

e(n)(f){n) — — / dm (g D + ih D sign(n — m))0(m), (34) 

2 Ji 

where e(n),gD and /ip are given by eqs.()3ip. As shown above, eq.fHJ) is satisfied by an 
eigenstate of the RD model with Ere = 0, and can be solved in the same way as was done 



in references 



22| . Taking the derivative with respect to n yields 



^ [e(n)0(n)] = ih D 4>( n ), (35) 



whose integral determines the functional form of 4>(n) 



A / f n dn' 
n) = — — exp ih D / — — - . (36) 



£{n) \ J 1 £{n' 
Using eqs.(|3*T|) and (J52|) this leads to 

, / % A (%E [ n dn' \ 

which coincides with the wave function for the generalized BK model. This result 

agrees with the fact that the commutator [X, P] converges asymptotically to ih (see figEJ). 
To find the eigenenergies E — 1/Ei let us consider eq.([34|) at the boundaries of the interval, 
n = l,N 



i r N 

£(1)0(1) = -(g D - ih D ) I dn<j>{n) (38) 

TV 



e(N)<j){N) = -{g D + ih D ) / d^(n) 



Dividing both eqs, 



9d + ihp = e(N)(P{N) 

g D -ih D £(1)0(1) ' 1 } 
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and using one finds 



9D + ih D f f N dn , 

exp t/ij, / — - , (40) 



or equivalently 



gj — ihj \ h J x e(n 

The eigenenergies are obtained taking the log 



9i + ihi (iE f N dn . 

ex P Hr / • ( 41 ) 



. . E f dn ot , , 

Ni(E) = — - / — - - -, 42 
v 7 2txH J 1 e(n) tt' v ; 



a = Arctan ( — ) (43) 



with 



and where Nj = 0, ±1, ±2, . . . label the eigenstates. In the case e(n) = n,gj = 0, (|42|) 
reduces to 

Ar,( £ ) = JL,ogAr-I (44) 
which yields an equally space spectrum symmetric around zero energy, 

ihi|> +1/2) - --o.* 1 .- («) 

In the limit N oo, this spectrum becomes a continuum with constant energy level density 
The conclusion is that the regularized BK model does not have a discrete spectrum whose 
eigenenergies could be identified with the Riemann zeros or an approximation to them. On 
the other hand, eq. (|44jl agrees with the leading term in Connes's formula (jlljl . which also 
diverges with the cutoff A. 

The wave function (|37)l in the case e(n) = n is 

^(n) = nl/2 _ iE/h , (46) 
which agrees with the BK wave function (|13|) . This is the wave function of a free particle 
moving in a box of length = log A. Indeed, let us define the coordinate 

q = log n, < q < L N = log A, (47) 
15 



and the momentum variable 



E 27T 1 

k = j, K = T^ i ~ n+ 2 ) - (48) 

The scalar product of two eigenstates (j4*E|) with quantum numbers n\ and ri2 is 

f N dn C Ln 
{^ Eni \iPeJ = A 2 / -ULn^-E^yn = A 2 dq = s (4Q) 

1 Ji n Jo 

— 1/2 

where A = L N is a normalization constant. Notice that the plane waves satisfy antiperi- 
odic boundary conditions, which can be changed by choosing gj ^ 0. In particular, gi — oo 
yields periodic BC's. 

The variable q can be defined for generic choices of levels e(n) 

r n dn' f N dn' 

?(n)= / 0< g <L N = / -£L (50) 

The momenta k is still defined by (148)1 . The eigenenergies £7, or rather momenta k = E/h, 
satisfy (jUJ), which can be written as 



e 2ta = e tkL N _ ^ 

This is the quantization condition of a free particle moving in a box of length L^, with 
twisted boundary conditions fixed by a. We give below a RG interpretation of q. 



IV. RENORMALIZATION GROUP ANALYSIS 



The Hamiltonian Hj is renormalizable in the sense that upon integration of the high 
energy degrees of freedom, the effective Hamiltonian governing the dynamics of the remaining 
variables coincides with the original one parameterized by new coupling constants. For a 
discrete QM model the RG procedure consists in the Gauss elimination of the highest energy 
degree of freedom and its replacement into the Schrodinger eq. for the other variables (one 
can also eliminate the lowest energy component in which case the flow goes towards the 
ultraviolet) |20j. Let us describe this process for Hj. We start from eq.()27j) written as 

(El - -gifl)lp n = fnfm^Pm + " fnfm^ra- (52) 

m=l m=n+l 
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Next we eliminate tp^ in terms of ifj n 



<N 



N-1 



TpN = Jj + lkl f2 52 fNfm^Pm- (53) 

2Ej - gi f N ^ 

Plugging this equation back into (|5^j). for n < N, one obtains a system of equations for the 
variables ip n <N, which is identical to the original system (jo^j) except that the couplings are 



a 2 + h 2 

g< gi+ Si + lh (54) 

hj = h I} (55) 
f n = f n (n = l,...,JV-l). (56) 

Hence /i/ and / n are RG invariant couplings, while gj changes under the RG. This is the 
reason to add the gj coupling to the Hamiltonian Hi since it is generated by the RG. 



The RD model ()28|) leads to an equation similar to eq. ()54|) for g^, where the term Ejf 



-2 



N — 

EiEn is replaced by En — -Erd, with E^j) the energy of the eigenstate. To study the low 
energy degrees of freedom, i.e. \E RD \ « e n , one can replace e^ — E RD by En, obtaining 
a RG equation for g^ which does not depend on -Erd- However, to analyze eq . (|54j) in 
the continuum, we cannot eliminate the term Ej, otherwise the scale dependence would 
dissapear. We shall use instead the parametrization (pTTj) to express gi,hj, fjsr in terms of 
the couplings go, h D ,EN of the associated RD model. Eq. (p>4"|) implies for g D 

J „ i 9d + h 2 D , , 

9d = 9d + , (57) 

zen — gD 

where go and hp must be regarded as functions of Ei, or E = 1/Ei. Eq. (JH7|) is the 
exact RG equation associated to a state with zero energy in the RD model. This result 
is a consequence of the relation between the X and RD models explained in section III. A 
pictorial representation of eq. ()57|) is given in fig. HJ 

Let's analyze eq. ()57|) in the continuum limit as in the previous section. For weak couplings 
we get 

g D (N-l)^g D (N) + 9 -M^l±^, (58) 

ZE N 
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FIG. 4: Graphic representation of ea. (|57j) . The term g 2 ^ + h 2 D arises from the tunneling of a pair 
of electrons from a level n.j < N to the highest level iV, followed by the decay to a level rif < N. 
The denominator 2eAr — 9d is a kinematical factor. 

where = go{.N) and ^ = go(N — 1) are regarded as a continuous function of N. A 
gradient expansion of go(N) yields the differential eq. 



(59) 



tfcg _ 9 2 p + h 2 D 
dN 2e N 

which for the uniform model becomes 

t = &° +h ^ (»> 

where iV(s) = e~ s N is the system size at the scale s and iV(0) = N is the initial size of the 
system. In general, for a monotonically increasing function e(n), one can define the scaling 
variable s as 



e(n') 

so that eq.(|59*|) takes also the form of eq. (jrjUj) . Comparison of eqs.floT"]) and (j5Uj) yields 



S W = / TTCTv (61) 



s(n) = L N - q(n), (62) 
which relates s to the variable q used in the previous section. The solution of (|6U|) is 

9d{s) = h D tan (\h D s + Arctan {^J~^j ^ > (63) 
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therefore 



9i{s) = i tan s + a^j . (64) 
The coupling gi(s) is periodic under the RG with a period 

2tt 2irh 

x ^h~ D -^r- (65) 

In the RD model the RG period of all the flows is the constant 2ir/hD- However, Ijfifij) . yields 
a RG period which depends on the energy of the state. This fact is related to the energy 
spectrum as we explain below. 

Russian doll scaling: gapped versus gapless 

Let us consider the n th state E n (N) given by eq.( 



E n (N) f N dm a 
2nn J x £[m) n 

After a RG cycle the size is reduced from iV to iV(A n ) (recall eqs.(|B3|l and (|6T])). 

2nh r N dm 



K E n (N) J N[Xn) e(mY (6?) 



so 



^ JN(\ n ) e{m) 

Splitting the interval of integration (1,-AT) in eq. ()66j) into the intervals (l,iV(A n )) U 
(N(X n ), N), and using (|55j) one gets 

EJN) f N{Xn) dm a 
2vrn A e(m) vr v 7 

This equation has the same form as eq. ()66|) with the replacements n — > n—1 and iV — > N(X n ), 

which implies, 

£ n (iV) = ^_!(Ar(A n )). (70) 
In the uniform case e n = n, gj = one has 

lo£f N „ n-1/2 

An = ^+I72' AT(A n ) = e- A "iV = iV™, (71) 
19 



(5a) 



(5b) 



FIG. 5: 5a) Energy levels of the RD model. The bound states lie below the lowest energy e\, while 
the antibound states lie above the highest energy e^. There are also states within the energy band 
(ei,£jv) ( n °t depicted). The scaling is exponential (see ea.([74j)). 5b) Energy levels E^ n = 1/E n of 
the I model (see ea.(|45|)). The energies converge to the energy e± = as ~ from above and 

below. 

so 

n-l/2 

E n (N) = E n _ 1 (N*V3), (72) 

which can be readily verified using eq. (|45j) . In the usual RD model the RG period is the 
constant = 2n/hD, which does not dependent on the energy of the states. The analogue 
of scaling relation (|72jl is 

E n (N)=E n - 1 {e- XD N). (73) 
This relation leads to an exponential decaying behaviour of the bound state energies 

E n (N) ~Ne- nX °, (74) 

which is in sharp contrast with the linear decaying of the energies Ej n = 1/E n in the 
X-model (see eq. (JH3)), which reflects the gapless nature of its spectrum (see figEJ)- 
Counting RG cycles 
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In a finite system the number of RG cycles is of course finite. In the RD model this 



number is given approximately by |21| 
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n c ~^logiV. (75) 

Z7T 

The reason for (|75jl is that after each cycle the system size is reduced by a factor e~ x ° 
where Xd = 2n/h D . Hence, the number of cycles n c needed to reduce the system to one 
site satisfies e~ ncXD N ~ 1, which leads to eq. (|75j) . This eq. also gives the number of bound 
states of the model, due to the one-to-one correspondence between RG cycles and bound 
states (recall (J73J)). 

For the uniform X-model we expect the formula (J75)) to give the number of RG cycles for 
the appropiate value of ho as a function of E n , namely 

n c (E n ) ~p!-\ og N = n + 1/2. (76) 
Iterating eq.(J?2J) one can compute the size of the system after n c RG cycles, 

n— 1/2 n-3/2 n+l/2-n e 

j\T jVn + 1/2 _> jYn + 1/2 n + l/2 \ (JJ) 

Hence the integer n, labelling the state E n , coincides roughly with the number of RG cycles 
n c needed to reduce the system to one site. This result is exact to leading order in N, but 
as we shall see below there are finite size corrections which will play an important in the 
discussion. 



V. EXACT SOLUTION 

The discrete Schrodinger eq.(|27|) can be solved in a very simple way which parallels 
the derivation done in the section III. This solution coincides with the exact Bethe ansatz 



solution of the RD model found in reference 



for a single Cooper pair. 



here is another 
The idea is 



derivation of the exact solution using the RG method of Glazek and Wilson 
to keep the energy E in the Gauss elimination procedure until the system size is one. That 
gives an exact equation for the eigenenergies E. The renormalizability of the model makes 
this procedure doable and that is the reason for its solvability. 
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As explained in section II, each state with energy Ei of the Hamiltonian ()19|) coincides 
with the zero energy state of the RD Hamiltonian ()25|) under the identifications (p?Tj) . The 
equation of that state is given by with E RD = 0, i.e. 



N 

1 - 

e 



""" ~~ \ ^2( gD + ih ° si S n ( ra ~ (78) 



2 

m=l 

Substracting the eqs. for (p n and (p n +i 

e n+1 (j) n+1 - e n (f) n = -rp(0n + (79) 

gives a discrete differential equation, which is a recursion relation for as a function of 

e n ,e n+1 and n , 

£ « + ^/ 2 n = l,...,iV-l. (80) 



4> n e n+1 -ih D /2' 
A gradient expansion of n and e n in (f?9*j) reproduces eq.(|3l)j). Iterating (jHUJ) yields 



TT e n + ih D /2 



J1 „ 



, e^-ihn/2- (81) 



Eqs. f|8U|) is a set of iV — 1 equations while ((75)) contain iV equations. Hence there is one 
more equation to impose. It is convenient to choose eqs. ()78|) at the two boundaries, n = 1 
and N obtaining, 



1 N 

(ex - ih D /2)(f) 1 = -(g D - ih D ) (82) 



m=l 

N 



Dividing both eqs. 



{e N + ih D /2)<p N = -(g D + ih D ) ^ m . 

m=l 



9d + ihp _ {sn + ihp/2) <f) N 
g D - ih D (ex - ih D /2) 0i 



(83) 



and using (jHH) gives 



9d + ihp _ jt- e n + ih D /2 ^ 

n=l ' 
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which can also be written as 



gi +ihi _yr £n + iE/2h 

gi - ihj ~ AA e n - iE/2ti [ ' 

where E is the energy of the BK Hamiltonian. The solutions of eq. ()85|) are obtained taking 
the log 



7T 



where Nj is an integer labelling the eigenstates. 
1) Completeness of the spectrum. 
Let's write eq.(f8*oT) as 



E \ a 



1 / 
Ni(E) = -J^Arctan I 

71 n=l \ 



(87) 



2he„ ) 7r 

n=l 

The condition e n = f~ 2 > 0,Vra implies that Nj(E) is a monotonically increasing function 
of E varying in the interval 

N a N a 

~ ~ ~ < ME) < y - J (88) 

where — oo < < oo. There are exactly iV integers in the interval (|8*Hj) corresponding to 
all the eigenvalues of H 1 ^. For some choices of N and a one may eventually find a solution 
with = oo. If iV is even all the solutions are finite. 
2) The spectrum of the operator P 

A particular example of BK Hamiltonian is given by the choice v n = e n = 1, gi = 0, which 
corresponds to the operator P (see eq.()22jl) . Its exact eigenvalues, p n , follow inmediately 
from eq. (|8*7j) . 



p n = 2 tan 



7T 1 



N N 

J < n < -, (89) 



which in the continuum limit coincide with the momenta of a free particle in a box with 
antiperiodic BC's, p n = j^(n + 1/2). 

3) Relation with the continuum approximation. 
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The power expansion of the function Arctan in ()87|) is, 

N W = £nti-- w +0(E>). (90) 

n=l 



In the continuum limit 



'» dn 



£n Ji e(n) ' 



(91) 



n=l 

one recovers eg. (|42j) . The higher order powers of E in (18 7j) are multiplied by terms of the 
form X]l/ e n +2m with m > 1. In the case where e n = n only the sum ^2,1/ e n ~ log AT 
diverges with AT, which controls the large N limit of Ni(E). 
4) Continuum approximation of the exact solution 

The previous discussion suggests to take the continuum limit directly in the exact equation 
(JS7J), namely 

1 f N ( E \ a 
NAE) ~ - / dn Arctan - - . (92) 

In the uniform case, e(n) = n,a = ir/2, and in the limit N >> \E\/2h » 1, eq. ([9*2*)l 
becomes 

^'li^-siW- 1 )^ 1 '' (93) 

The finite piece of this equation agrees with Connes's formula (JUJ, except for a term linear 
in namely E/(2tv) logrr. The origin of this term, in the Riemann counting formula, can be 
traced back from the expression (|7J) . and it is due to a factor vr 5 appearing in the functional 
relation satisfied by C( s )- The leading term E/(2tv) logE, as well as another linear term in 
E, come from the log of the Gamma function in (JJJ). 

The divergent terms in eqs. ()93|) and (fTTj) have a similar form which depend on the re- 
spective cutoffs iV and A. Let us suppose for a while that both cutoffs are related by the 
eq. iV = 7rA 2 . Then the divergent and the finite parts in both formulas agree. Notice 
that the n factor in the latter relation explains the missing factor E/(2n) logn in eq. (|93jl . 
Unfortunately, the relation N = 7rA 2 does not seem to follow from the counting of states 
in both models. In the X model this is given by N, while in the Connes model, at the 
semiclassical level, it would be given by A 2 /7r (h = 1) corresponding to the phase space 
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(0 < x < A, —A < p < A). This comparison gives A" = A 2 /tt, rather than A" = 7rA 2 . In 
any case, the linear term oc E, in both counting formulas, depends on the cutoffs and hence 
on the particular regularization choosen. Apparently, Connes's regularization and ours are 
different. In the next paragraph with shall consider a zeta function regularization of the 
model, which will shed further light on this issue. 

Another interesting point concerns the spectroscopic interpretation of our results. As we 
said above they are along the lines of Connes's absortion picture. However it must be kept in 
mind that the eigenstates counted by the smooth part of the Riemann formula are not really 
missing in the whole spectra but shifted to higher energies because the interactions. This 
blueshift makes that in a range of energies, say (0, E), there are less states than expected 
from the analysis of the continuum limit. In this sense, a more appropiate spectroscopic 
interpretation of our results will be in terms of a blueshift of energy levels, which are then 
missing in fixed energy intervals. 

5) Exact solution in the uniform case. 

Eq.fjHZj) can be given an exact analytic formula in terms of known functions in the uniform 
case. We shall add a zero point contribution, a, to the energy levels, i.e. 

e n = n + a, n = l,2,...,N, (94) 

which does not modify the large N properties of Nj discussed above. It is more convenient 
to write (|87jl in the product form (J85|) (h = 1) 

N 

2mNi(E) = -2ia TT gn + iE/2 

L\e n -iE2 1 ; 

rc=l ' 

Inspired by (J93|) . we shall define the finite part of Nj as 

ni(E) = lim (JLlogN-NjiE)), (96) 

iV— >oo \Z7l J 

where a minus sign has been introduced to take care of the relative minus in eq.(|93|). The 
expression for nj(E) follows from (|93j) . 



N 

e ~iE log N e ~2ia 

Using the eqs. 



N 

-2 mni (E) = Um -iElogN -2ia TT £ » + iE / 2 (Q7) 
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I 00 

n(i+s)«-^. («) 

v ' n=l 



7 

AT 



one finds 

N 

e -iE log N 

N—*0O 

which plugged in (|57jl. and taking the log, gives 



lim e -, W TT n + q + ifl/2 = r(l + a-^/2) 

£A n + a - r (1 + a + iE/2) ' v 7 



1 / £\ a 

n 7 (E) = -logr l + a + i- + -. (101) 
7r y I J n 

Comparing this equation with the smooth part of the Riemann's counting formula (J7J), we 
deduce that the Gamma term in both eqs. agree if we choose 

a = -| (102) 

The energy levels e n (see eq. (I94J) ) corresponding to this choice are those of a harmonic 
oscillator with zero point energy 1/4, instead of 1/2. This vacuum energy arises from 
Neumann BC's at the origin which select the even eigenf unctions under parity. In this sense 
the RD model can also be seen as a harmonic oscillator with Neumann BC's perturbed by 
the gr> and hn interactions given in eq. (|25|) . 

The asymptotic expansion of (jlOlj) gives the finite part of (jHHJ), and it is related to the 
smooth part of the Riemann formula (J7J), for a = 7r/2, as 

N sm (E) = ni {E)-^\og-K + \. (103) 

Z7T Z 

Notice again the term — i?/27rlog7r which, as we said, is regularization dependent. To 
highlight this point, we shall apply a zeta function regularization to the expression (|97j) . 
Choosing N = 00 one has formally, 

-2mni(E) _ -2ia Tln=l M £ n + iE/2) 
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where —ni(E) stands for the regularized value of Ni(E) and p, is a regulator. The infinite 
products in (jlU4|) can be regularized using the zeta function (37 1 

Hv(n + z)=^ (-j=r(z)) . (105) 

n — n \ v / 



The result is 

e -2mn l{ E) = & -2ia -iE T(l + a-zE/2) 

P r (1 + a + iE/2) ' v ; 

and 

1 / E\ ot E 

n I {E) = -\ogY (l + a + i-) + - + — log/i. (107) 

7T \ Z y 7T Z7T 

Choosing \i = l/n one gets the factor — _E/27rlog7r, and the relation ()103|) is replaced by 

N sm (E) = ME) + \- (108) 

There is also a mismatch of 1/2 due to the a term, which is not important for large values of 
E but which can be relevant for small ones. We have carried out a numerical computation to 
assess the accuracy of N sm (E n ) for predicting the position of the zeros (E n is the imaginary 
part of the n th Riemann zero). The closest N sm (E n ) comes to n, the better the approximation 
is. The results for the first 10 zeros are collected in table 1. They show that N sm (E) + 1/2 = 
nj(E) + 1 is a much better fit than N sm (E), suggesting that the zeros are associated to 
complete RG cycles. This result is confirmed in figure |H1 where we plot the difference 
n — N Bm (E n ) — 1/2 for the first 40 zeros. This extra factor 1/2, that improves the location 
of the zeros, was also obtained by Berry j(| and Badhuri et al. j^cj in their respective 
approaches. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


E n 


14.1347 


21.0220 


25.0109 


30.4248 


32.9350 


37.5861 


40.9187 


43.3270 


48.0051 


49.7738 


N sm (E n ) 


0.4497 


1.5702 


2.3936 


3.6710 


4.3172 


5.5935 


6.5651 


7.2943 


8.7708 


9.3483 


N sm {E n ) + \ 


0.9497 


2.0702 


2.8936 


4.1710 


4.8172 


6.0935 


7.0651 


7.7943 


9.2708 


9.8483 



Table 1.- Values of the smooth part of the Riemann counting formula, N sm (E n ) for the first 
10 Riemann zeros E n (eq.(J7j)). N sm (E n ) + 1/2 is also given for comparison. 
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5 10 15 20 25 30 35 40 

Order of the Riemann zeros 



FIG. 6: Deviation of the first 40 zeros. The mean and square mean of n — N sm — 1/2 are given by 
-0.003825 and 0.188102 respectively. 



The oscillations in n — N sm (E n ) — 1/2 are randomly distributed due to its relation to 
C(l/2 + iE). This fact implies that the uniform X-model, eventhough describes the smooth 
part of the Riemann's zeros, does not explain the origin of their randomness. 

To solve this problem we have tried to modify the uniform energy levels ()94|) in several 
ways. The first one is to modify the position of the energy levels. A slight modification of 
M|) consist in the addition of a 1/n correction, 



£ n = n + a + -, 
n 



n=l,2,...,iV, 



(109) 



which does not change the logiV behaviour of Nj. Using (}9~Tj) one readily finds 



m(E) = - log [r (i + a+ (E)) r (i + o_(£))] + -, 

7T 7T 



(110) 



where 



a±(E) 



1 
2 




a + i- 



E 



4c 



111^ 



For large values of E, eq. ljllOj) agrees with eq. (j9Hjl to order ElogE, but the comparison 
with (J7J) is lost since we have now the product of two gamma functions instead of one, as in 
(llOlj) . Adding higher order corrections to ()109j) of the form 1/n 2 , . . . gives similar results. 
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Another possibility is to eliminate some energy levels in e n , for example those associated 
to the prime numbers. This choice is suggested by the quantum chaos hypothesis accor ding 
to which the primitive orbits of the chaotic Hamiltonian are labelled by the primes 



14. 



15|. The truncation of the prime energy levels, i.e. e p (with p a prime), changes the 



asymptotic expansion of Nj which behaves as log(iV/ log N)) instead of logiV. However its 
finite part, nj, does not improve the location of the Riemann zeros. These negative results 
can be understood from the formula (|57jl for Nj. Indeed, the function Arctan(.E/2e n ) 
varies smoothly between and tx/2 in a range of E set by e n . The same applies for a 
superposition of those terms with different energies, making very difficult to obtain a random 
curve interpolating the zeros. The uniform energy levels seems to provide the best possible 
approximation within the X model. This fact has lead us to a further generalization of this 
model. 



VI. THE T± MODELS 

The couplings gi t D have so far played an auxiliary role in the construction of the X 
model. They appear in the renormalization and the exact solution, but they did not take 
part in the dynamics except for setting the boundary conditions. This suggests that giD 
must play a more significant role in the dynamics underlying the Riemann zeros. A natural 
generalization is to replace gj, in the Hamiltonian (jl9|) . by a generic real symmetric matrix 
gnmi however the hj interaction, no longer keeps its simple form under the RG, becoming a 
generic matrix. Surprisingly enough, there are two choices of g l n m , which leave the hj term 
invariant under the RG. The corresponding Hamiltonians are given by, 

(n\H I± \m) = -f n (g I± {n,m) + ilu sign(n - m)) f m , (112) 

where the matrix elements gi ± (n, m) are defined in terms of a set of couplings gi, p (p = 
1, . . . , N) as 

{max(n, m) for X + 
(113) 
min(n, m) for X_ 

The structure of gj ± (n,m) is displayed in fig. [7[ which also shows the direction in which 
the models are renormalized. For iV finite the X± models are related by the transformation 
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FIG. 7: Graphical representation of the coupling matrices gi±(n, m) given in ea. ()113j) . The 
continuous lines represent common values of g p = gj tP . The labels n, m (and the energies e n ) 
increase from left to right and from top to bottom. The RG in the X + model runs towards the 
infrared (IR) and in the X_ model towards the ultraviolet (UV). 

n —> N — n and hj — > —hi, which changes the order of the energy levels e n . However in the 
limit V — > oo they are inequivalent. We shall assume, as usual, that e n increases with n. 

The correspondence between eigenstates of the Z and the RD models is also maintained 
for Z±, denning the associated RD ± models as 

(n\H RD± \m) = e n 5 ntm - - (g D± (n,m) +ih D sign(n - m)) , (114) 

where gn^n^m) has the form (|113|) . with g Ip replaced by gr>,p- The relation is established 
between an eigenstate, with energy Ej, of the Z± models and the zero energy state, Erd± = 
of the RD ± , with couplings constants given by eqs. (}3Tj) . where the equation g D = gi/Ej is 
replaced by gr,, P = gi, P /Ej. If gj, p = gj, \/p the I± models reduce to the original Z model 
and the same applies to the RD± models which become the RD one. 

From a physical viewpoint the Z + and RD + models are characterized by the fact that the 
tunneling among two energy levels, induced by the g couplings, depends only on the state 
with the highest level n, i.e. the highest energy one for RD + . For the Z_ and RD_ models, 
it is the lowest level that matters. Looking at fig 0] one realizes that these models have to be 
renormalized in different ways for h^o to be invariant. In the RD + model one has to eliminate 
the highest energy mode, n = N, which involves the product {gD,N J rih£,){g£, t N — ih£ ) ) which 
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is real, leaving hp invariant. However, in the RD_ model, the level to be eliminated is n = 1, 
which also lives ho fixed. Hence the RG flow goes towards the infrared in the RD + model, 
and towards the ultraviolet for RD„. Of course the RD model can be renormalized in both 
ways. We shall not write explicitely the RG eqs. for the couplings but they can be easily 
derived using the techniques of section IV. 

A. Continuum limit 

Let's first consider the RD + model which, as we said, corresponds to T + . The Schodinger 
eq. (jll4|) . for a state with Erd = 0, is 

^ n N N 

£ n <Pn = r ^ gp,n4>m + 7, gp,m4>m + ihp \, sign(n - m)(f) m , (115) 

m=l m=n+l m=l 

and in the continuum 



1 f n 1 f N ifi£> f N 

e(n)(j)(n) = - / dm gpin) 0(m) H — / dm ^r)(m) 0(m) H / dm sign(n — m) <p{m). 

2 Ji 2 J n 2 J 1 

(116) 

The derivative respect to n 

d 1 dgD f n 

— [e{n)4>{n)\ = ih D <p{n) + J dm<p(m), (117) 

is an integro-differential equation for <f)(n), subject to the boundary condition 

1 ^ 



e(N)<j>(N) = -(g D (N) + ih D ) J dm <f>(m). (118) 
Eq. (jll7J) can be converted into a second order differential equation for the function 

x( n ) = / dm 4>(m), (pM = — , (119) 
Ji dn 

namely 

d (c(n) dx ^\ ih u dx ld9D x = (120) 
dn V dn / dn 2 dn 
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together with the BC's 



X(l) = 0, e^(N) = \{g D {N) + ih D ) x (N). (121) 
Making the change of variables n — > q (recall eq. (|5()j) ). the eqs. (|12()jl and ()121j) become 

d 2 X .. dx ldg D 

d?- %kD dj-2lj- X = > (122) 
dx 1 

X = (at q = 0), — = -(g D + ih D )x (at q = L N ). (123) 

dq 2 



The second summand in eq. ()122j) can be eliminated by the gauge transformation, 

X = e~ ihDq/2 x, (124) 



leading to 



d 2 x 1 dgn ~ h 



X + ^X = 0, (125) 



dq 2 2 dq ^ 4 
dy 1 

X = (at q = 0), = x (at g = (126) 

Notice that ()125j) is the Schrodinger equation of an effective Hamiltonian 

corresponding to a potential given by the derivative of the function goio)- The energy of 
that state is given by E + = h 2 D /4. It is interesting to rederive the results obtained in previous 
sections where go is constant, i.e. V+ = 0. The solution of eq. ()125j) gives a superposition of 
plane waves 

X = e ihm/2 + Ce' lhDq/2 , (128) 

where C = — 1 to satisfy the BC at q = (eq. ([126|) ). The BC at q — gives the plane 
wave quantization (see eqs. pOj) and (fHT|l ) 

e i hD L N = 9D + ihD (129) 
g D - ih D 
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The RD_ model can be studied in a similar manner. We give for completeness the results. 
The function x( n ) m ust be defined as 

f N d X 
xi n ) — / dm 4>(m), <f)(n) = (130) 

X is still given by eq. (jl24j) . The effective Hamiltonian for x is 

+ (131) 
and the eigenenergies are given by the same formula E_ = h 2 /A. The BC's also change 

dx 1 ~ 

= —j9\ X (at q = 0), x = (at q = L N ). (132) 

As in the X model these results can be generalized to the discrete case, obtaining an exact 
equation for the eigenvalues of the 1± hamiltonians. The analogue of eq.(|85p is a matrix 
like Bethe equation, related to the fact that the wave function Xn satisfies a second order 
discrete differential equation. The results will be presented elsewhere. 

Eqs. (I122|) and (|131|) bring the idea of a potential V + or V"_, whose scattering theory would 
produce the oscillating part of the Riemann counting formula. Pavlov and Fadeev showed 
long time ago that the zeta function ((s), on the line Re(s) = 1, appears in the scattering 



phase shift of particles moving on surfaces of constant negative curvature |38l. I39t 40]. This 
is a possible direction of research, which is likely to be related to quantum chaos. On the 
other hand, our approach is based on discrete Hamiltonians, so the solution of the previous 
problem is a possible strategy to find the correct choice of the discrete couplings g n leading 
to the precise location of the Riemann zeros. 

On more general grounds, let us recall that one of the motivations to consider the Hamil- 
tonian xp was the breaking of time reversal symmetry, which should be related to the GUE 
statistics of the zeros. This Hamiltonian gives an accurate semiclassical description of them 
but not of their fluctuations. The I± models break time reversal, but they also break the 
reversal of the RG time direction. Could this additional breaking be related to the GUE 
statistics of the zeros? 

There are another interesting questions regarding the integrability of the models consid- 
ered in this paper. The RD model, including its many body version, is exactly solvable a la 



Bethe and integrable (i.e. infinite number of conserved quantities) 



3- 



Thanks to the map 
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X — > RD, we have been able to solve the X model and in turn the regularized BK model. 
Is the X model, including its many body version, integrable? It is not obvious how to gen- 
eralize the map X — > RD to the many body case, due to the existence of several rapidity 
variables. The QM T± models are renormalizable and exactly solvable, although the ana- 
logue of the Bethe equation has a more complex structure. In the RD model the coupling 
parameterizes a boundary operator in the transfer matrix of an inhomogenous vertex model. 
The existence of several couplings go, n suggests that the corresponding transfer matrix, if 
it exists, involves more than one boundary operator. 

VII. CONCLUSIONS AND PROSPECTS 

We have shown in this paper that the Berry and Keating Hamiltonian can be quantized 
in a consistent way on a lattice. This quantization has been achieved thanks to the relation 
of the BK Hamiltonian to two QM models which have a cyclic RG, specially the RD model 
of superconductivity. The later models are renormalizable and have an exact solution a la 
Bethe, which permits a detailed study of their spectrum. 

The first of these models, i.e. the X model, is the inverse of the BK Hamiltonian plus 
an extra coupling which, in the BCS framework, is the standard pairing interaction. This 
model has a continuum spectrum related to RG cycles with an energy dependent period. 
The second model is the RD model, which has bound states related by a scale factor which 
depends on a fixed RG period. These two models are intimately related, since an eigenstate 
with energy E of the X model, can be mapped into a bound state at the threshold of an RD 
model whose RG period is given by E. The map, which is exact, establish an unexpected 
correspondence between the BK model and the RD model of superconductivity where the 
time reversal symmetry is broken explicitely. 

Using these QM models we have given an spectral interpretation of the smooth part of 
the Riemann counting formula of the non trivial zeros. In the X model it counts the missing 
states in the continuum spectrum below a given energy. The result depends on the cutoff, 
which in our case is given by the number of sites. This result seems to agree with Connes's 
absortion spectral interpretation in the adelic theory, but the cutoff is different and a precise 
comparison with ours is not conclusive. In the RD model, the smooth part of the Riemann 
formula gives the number of missing bound states with respect to the leading term which 
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follows from the scaling properties of the cyclic RG. It is a finite size correction of the 
Russian doll scaling. In a certain sense, it can also be seen as an anomaly for the discrete 
RG transformations, i.e. the RG cycles, which leave the Hamiltonians invariant. 

We have looked for a choice of parameters of the X model that would explain the oscillating 
part of the Riemann formula. In the RD language this means choosing the energy levels. 
However a numerical study shows that the best choice is given by equally space energy levels. 
It thus seems that the origin of the random position of the zeros lies beyond the X model. To 
explain this randomness we propose a natural generalization of the X and RD models where 
the coupling g is replaced by a set of discrete couplings g n which depend on the level n. The 
guiding principle is renormalizability, which in this context means the RG invariance of the 
time reversal breaking interaction. Quite surprisingly, there are two models, X±, satisfying 
this condition. They differ in the way the RG procedure is implemented. In the X + model, 
the RG runs from the UV towards the IR, while in the X_ model the order is reversed. 

Finally, we have begun to explore the properties of the X± models in the continuum limit, 
finding two QM models with different potentials related to the gradient of the coupling 
function g n and different boundary conditions. We suggest that an appropiate choice of 
these potentials, and in turn of the couplings constants g n , may explain the local fluctuations 
of the Riemann zeros. The 1± models are likely to be related to other approaches to the 
Riemann zeros, specially to quantum chaos. It would be interesting to investigate that 
connection. 
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